AO- A 005  319  ILLINOIS  UNIV  AT  URB ANA-CHAMPAIGN  APPLIED  COMPUTATION— ETC  F/G  12/1 

PLANE-SWEEP  ALGORITHMS  FOR  INTERSECTING  GEOMETRIC  FIGURES* CU> 

OCT  79  J  NIEVER6ELT  *  F  P  PREPARATA  N00014-79-C-0424 

UNCLASSIFIED  ACT-10  NL 

ur  ■■■■■■■■■■■ 


9P_ 


ADA08531  9 


™  ■■ 1  a  - «a  aaiaaaauaaa 

COORDINATED  SCIENCE  LABORATORY 

APPLIED  COMPUTATION  THEORY  GROUP 


PLANE-SWEEP  ALGORITHMS  FOR 
INTERSECTING  GEOMETRIC  FIGURES 

J.  NIEVERGELT 
ER  PREPARATA 


D1W 

J^ELECTEgSL 


m  JUN  9  1980  •  9 


REPORT  R-863 


ose  and  «cd®' 


ase  anu  - 

ttriirtnWd* 


►  UNIVERSITY  OF  ILLINOIS  -  URBANA,  ILLINOIS 


.SECURITY  CLASSIFICATION  of  This  pace  ,'>Wi«n  Dot,  Enttrtd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


a.  recipient's  catalog  number 


4.  title  (And  SubtllU) 


PLANE -j^WEEP  ALGORITHMS  FOR  INTERSECTING 
(fEOMETRIC  FIGURES^  ^ 


7.  author?*; 


J.  Nievergelt  end  F.  P.  Preparata 


S.  TYPE  OF  REPORT  a  PERIOD  COVERED 


(i» ; 


ACT  - 1 8.U  ILU-ENG  78-2256 


a.  contract  or  grant  numberc*; 

NSF  MCS  78-13642; 
N0OO14-79-C-G424 


9-  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Coordinated  Science  Laboratory 
University  of  Illinois  at  Urbana -Champaign 
Urbana,  Illinois  61801 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Joint  Services  Electronics  Program; 
National  Science  Foundation 


14.  MONITORING  AGENCY  NAME  *  AOORESSflf  dllttrDt  /ran  Controlling  Office; 


l«.  DISTRIBUTION  STATEMENT  (ol  (fill  Report; 


T.  °°o  I  cr~r  T  ter 
AREA  A  WORK  UNIT  NUMBERS 


13.  NUMBER  of  paces 
22 


IS.  SECURITY  CLASS,  fof  rtle  report; 


UNCLASSIFIED 


I  So.  OECLASSI  FI  CATION/  DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  thm  ebstrmet  entered  In  Block  20,  If  different  from  Report) 


t9.  KEY  WORDS  (Continue  on  reverse  side  ll  necessmry  end  Identify  by  block  number) 


Computational  geometry,  concrete  complexity,  combinatorial  algorithms, 
intersection  problems,  maps,  polygons 


4  rtf  -  ’ 


20.  ABSTRACT  (Continue  on  reverse  side  If  necessmry  end  identify  by  block  number) 

We  considers various  types  of  geometric  figures  in  the  plane  defined 
by  points  connected  by  straight  line  segments,  such  as  polygons  (not 
necessarily  simple)  and  maps  (embedded  planar  graphs);  and  the  problem  of 
computing  the  intersection  of  such  figures  by  means  of  a  fgreedyf  type  of 
algorithm  that  sweeps  the  plane  unidirectionally .  Let  n  be  the  total 
number  of  points  of  all  the  figures  involved,  and  s  the  total  number  of 
intersections  of  line  segments.  .We  show-  that  in  the  general  case  (no 
convexity)  a  previously  known  algorithm  can  be  extended  to  compute  in  time  - 


DO  ,  1473  eomoN  OF  I  NOV  as  is  obsolete  UNCLASSIFIED 

,  •  y-..  /  SECURITY  CLASSIFICATION  OF  THIS  PAGE  flWipn  Dwa  Inter**; 

-■  U  O'-  / 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  this  FAOSOFXm  Dtm  Knft+4) 


PLANE-SWEEP  ALGORITHMS  FOR 
INTERSECTING  GEOMETRIC  FIGURES 

by 

J.  Nievergelt  and  F.  P.  Preparata 


This  work  was  supported  in  part  by  the  National  Science  Foundation 
under  Grant  MCS  78-13642  and  Joint  Services  Electronics  Program  under 
Contract  N00014-79-C-0424 . 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of 
the  United  States  Government. 


This  report  Is  issued  simultaneously  by  the  Coordinated  Science 
Laboratory,  University  of  Illinois  at  Urbana-Champaign,  and  by  the 
Ins ti tut  fur  Informatik,  Eidgenossische  Technische  Hochschule,  Zurich. 


Approved  for  public  release.  Distribution  unlimited. 


PLANE -SWEEP  ALGORITHMS  FOR  INTERSECTING  GEOMETRIC  FIGURES 


J.  Nievergelt  and  F.  P.  Preparata 


Abstract 

We  consider  various  types  of  geometric  figures  in  the  plane  defined 
by  points  connected  by  straight  line  segments,  such  as  polygons  (not 
necessarily  simple)  and  maps  (embedded  planar  graphs);  and  the  problem 
of  computing  the  intersection  of  such  figures  by  means  of  a  "greedy" 
type  of  algorithm  that  sweeps  the  plane  unidirectionally .  Let  n  be  the 

total  number  of  points  of  all  the  figures  involved,  and  s  the  total  number 

of  intersections  of  line  segments.  We  show  that  in  the  general  case  (no 
convexity)  a  previously  known  algorithm  can  be  extended  to  compute  in  time 

0  ( (tri-s)  logn)  and  space  0(n+s)  all  the  connected  regions  into  which  the 

plane  is  divided  by  intersecting  the  figures.  When  the  regions  of  each 
figure  are  convex,  the  same  can  be  achieved  in  time  O(nlogn+s)  and  space 
0  (n) . 
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1 .  INTRODUCTION 

One  type  of  algorithm  that  emerged  from  recent  research  in  computational 
plane  geometry  promises  to  be  efficient  for  all  problems  to  which  it  applies. 
It  sweeps  the  plane  "from  left  to  right",  advancing  a  "front"  or  "cross- 
section"  from  point  to  next  point.  All  processing  is  done  at  this  moving 
front,  without  any  backtracking,  with  a  horizon  or  look-ahead  of  only  one 
point.  Algorithms  of  this  type  are  often  called  "greedy,"  in  contrast  to 
exhaustive  search  algorithms  that  create  tentative  partial  results  perhaps 
to  be  discarded  later.  An  important  issue  in  computational  complexity  is 
to  understand  which  problems  can  be  solved  by  greedy  algorithms.  We 
contribute  to  this  issue,  and  present  two  efficient  algorithms  for  problems 
that  have  applications  in  geographic  data  processing. 

Shamos  and  Hoey  [6,7)  presented  an  algorithm  that,  by  sweeping  the 
plane  unidirectionally,  determines  in  time  O(nlogn)  whether  or  not  n  line 
segments  are  free  of  intersections.  Bentley  and  Ottmann  [2]  have  extended 
this  algorithm  to  report  all  s  intersections  of  n  line  segments  within 
time  O((s+n)logn) .  We  elaborate  on  this  type  of  algorithm  in  two  directions. 
First,  we  show  that,  within  the  same  asymptotic  effort,  it  can  compute  all 
connected  regions  into  which  the  plane  is  divided  by  the  line  segments, 
where  each  region  is  presented  by  a  cyclic  list  of  all  its  boundary 
vertices  and/or  segments.  Second,  we  show  that  assumptions  of  convexity 
allow  one  to  improve  these  asymptotic  bounds.  Specifically,  we  present  an 
algorithm  that  computes  all  s  intersections  of  two  convex  maps  (embedded 
planar  graphs  with  convex  regions)  with  a  total  of  n  points  in  space  0(n) 
and  time  O(nlogrrt-s),  which  is  optimal  in  s. 
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2.  REGIONS  FORMED  BY  A  POLYGON 

Later  sections  of  the  paper  will  refer  to  several  versions  of  the 

problem  of  computing  regions  of  the  plane  formed  by  embedded  graphs.  Here 

we  present  in  detail  the  version  which  is  simplest  for  expository  purposes: 

let  the  given  graph  be  a  polygon,  i.e.,  a  closed  chain  of  n  line  segments 

(or  equivalently,  a  cyclic  list  of  n  points  in  the  plane).  The  reader  who 

has  understood  how  the  algorithm  works  in  this  case  will  have  no  difficulty 

in  following  the  brief  description  of  how  it  applies  to  the  intersection 

of  more  general  plane-embedded  graphs . 

2.1  Statement  of  the  Problem  and  Terminology 

Given  a  sequence  of  n  points  =  (x^,y^),  i  =  l,2,...,n,  in  the  plane, 

a  polygon  with  vertices  is  the  sequence  of  line  segments  Vj^,  V2V3’'*'’ 

V^Vj.  These  n  line  segments  in  general  define  s  intersection  points 

Wj  =■  (Xj,y^>,  j  *  l,2,...,s.  When  s  =  0  the  polygon  is  called  simple  and 

divides  the  plane  into  two  regions,  an  internal  bounded  region  R^  and  an 
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external  unbounded  region  Rq .  In  general  s  =  0(n  ),  and  the  polygon 

divides  the  plane  into  r+1  ^  2  disjoint  regions,  namely  the  external 

unbounded  region  Rq  and  r  simply  connected  internal  regions  R^,...,Rr 

(when  the  polygon  is  non-degenerate,  r  =  s+1) .  Each  region  is  itself  a 

polygon  that  has  as  its  vertices  some  subset  of  [vi,...,V  ,W.,...,W  }. 

i  n  i  s 

The  desired  result  is  a  list  of  all  regions,  where  each  region  is  given 
by  a  cyclic  list  of  its  vertices,  starting  with  the  right-most  vertex;  the 
external  region  in  clockwise  order,  the  internal  regions  in  counterclockwise 
order.  Figure  1  illustrated  these  concepts. 


Fig.  1.  A  polygon  of  5  line  segments  and  the  regions  it  defines: 


R_:  V.W.V.V, 

0  4  115 

R.:  V.V.V^W^-V.W, 
1  4511321 


R2:  v2V3W1 


As  Fig.  1  suggests,  we  make  certain  assumptions  of  non-degeneracy, 
namely:  all  points  (originally  given  or  intersection)  are  distinct;  a  point 
lies  on  only  those  line  segments  on  which  it  must  lie,  and  no  others  (e.g., 
does  not  lie  on  V^V^) .  Section  4  discusses  the  modifications  required 
to  handle  degenerate  pictures. 

The  notions  above  are  sufficient  to  describe  the  problem  and  its 
solution.  In  order  to  describe  the  algorithm  which  computes  the  solution, 
we  Introduce  auxiliary  concepts  that  reflect  the  dynamic  aspects  -  a 
unidirectional  sweeping  of  the  plane.  Call  the  originally  given  points 
and  the  intersection  points  simply  points .  P^.P^, . . • sorted  in 
order  of  increasing  x-coordinate-  (We  assume  for  ease  of  exposition  that  no 
two  points  have  equal  x-coordinates ;  if  P  and  P.  .  have  equal  x-coordinate, 


a  lexicographic  ordering  on  the  pair  (x,y)  suffices  for  the  following 
discussion  to  apply.)  Under  our  assumptions  of  non-degeneracy  and  distinct 
x-coordinates ,  each  point  can  be  classified  uniquely  into  one  of  the 
following  4  categories : 

-  a  start  point  QCf  ' 

-  an  end  point  \ 

j  points  originally  given  as  vertices 

-  a  bend  J  of  the  polygon. 

-  an  intersection  poinOOC 


Fig.  2.  Figure  1  revisited.  Points  have  been  ordered  according  to 

increasing  x-coordinate .  Cross  section  y  contains  one  interval 
for  each  of  the  regions  R^  and  R2,  and  two  intervals  for  two 

branches  of  R^  that  have  split.  Oriented  "rubber  bands"  trailing 

the  cross  section  V  are  shown  as  broken  lines. 


A  cross  section  is  a  vertical  line  in  the  plane  along  with  all  the 


information  about  which  line  segments,  regions  and  components  it  cuts, 
and  in  what  order.  The  line  segments  cut  by  a  cross  section  partition 
it  into  intervals .  We  are  mainly  concerned  with  cross  sections  that  do 
not  pass  through  any  points;  the  set  of  all  (topologically  equivalent) 
cross  sections  that  lie  between  two  adjacent  points  is  called  a  slice . 

A  cross  section  that  passes  through  (exactly)  one  point  is  called  a 
transition . 

As  the  current  cross  section  sweeps  the  figure,  it  drags  along  "rubber 
bands"  that  hug  the  periphery  of  regions,  as  shown  in  Fig.  2.  Much  of  the 
specification  of  the  region-finding  algorithm  is  concerned  with  properly 
maintaining  these  rubber  bands  across  transitions. 

2.2  Data  Structures  Maintained  by  the  Algorithm 

The  region- finding  algorithm  to  be  presented  in  section  2.3  operates 
upon  three  data  structures.  These  are  crucial  for  understanding  the 
algorithm  as  well  as  for  its  efficiency.  We  present  these  data  structures 
at  a  fairly  abstract  level:  by  postulating  what  operations  must  be  performed 
and  how  much  time  is  available  for  them  (asymptotically  in  terms  of  the 
problem  parameters  n  and  s) .  We  refer  to  standard  textbooks  [1,4]  for 
concrete  implementations  that  realize  the  postulated  time  bounds. 

The  x-structure  X 

At  any  moment  X  contains  those  points  that  have  been  discovered  so  far 
and  are  yet  to  be  processed,  sorted  according  to  increasing  x-coordinate . 
Points  are  assigned  a  type  according  to  the  classification  of  section  2.1. 
The  x-structure  is  a  priority  queue  that  must  support  the  following 
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operations  within  time  bound  O(log  k)  when  it  contains  k  entries: 

-  MIN:  find  and  remove  the  entry  with  minimal  x-coordinate 

-  INSERT:  insert  a  new  entry  with  given  x-coordinate. 

Heaps  or  balanced  trees  are  suitable  for  implementing  priority  queues. 

Initial  content:  the  n  originally  given  points,  sorted,  with  their 
classification. 

Final  content:  empty. 

At  each  transition,  the  point  which  defines  this  transition  is  removed, 

and  at  most  two  intersection  points  are  inserted  into  X.  During  execution 

a  total  of  n+s  points  move  through  the  x-structure,  hence  the  maximal 

2 

number  of  entries  at  any  given  time  is  <  n+s.  Since  s  =  0(n  ),  any  operation 
on  the  x-structure  can  be  done  in  time  O(log(n+s))  =  O(logn). 

The  y- structure  Y 

Y  contains  all  the  information  about  a  cross  section  which  is 
representative  of  its  entire  slice.  It  has  an  entry  for  each  segment 
intersected  by  the  slice,  including  two  sentinels  corresponding  to 
y  =  +<»  and  y  =  -®,  and  thus  it  never  has  more  than  n+2  entries.  Y  is  a 
die tionary  (see  [1])  that  must  support  the  following  operation  within  time 
bound  O(logk)  when  it  contains  k  entries: 

-  FIND:  given  a  point  (x,y),  find  the  line  segment  in  the  cross  section 

at  x  whose  y-value  does  not  exceed  y. 

In  addition,  in  constant  time,  the  y-structure  supports  the  following 


operations: 
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-  DELETE(s):  given  a  pointer  to  an  entry  s,  delete  that  entry. 

-  INSERT (s ; t) : given  a  pointer  to  an  entry  s,  insert  a  new  entry  t 

following  that  entry. 

-  PREDECESSOR (s ): given  a  pointer  to  an  entry  s,  obtain  the  pointer  to 

the  immediately  preceding  entry. 

-  SUCCESSOR (s ): given  a  pointer  to  an  entry  s,  obtain  the  pointer  to 

the  immediately  following  entry. 

Various  types  of  balanced  trees,  with  additional  "prodecessor"  and 
"successor"  pointers  can  be  used  to  implement  such  dictionary.  Initial  and 
final  content  of  the  y-structure:  the  pair  (-■»,+*)  .  Figure  3  shows  the 
y-structure  of  the  slice  between  and  P^  in  Fig.  2. 


Pointer  -o  Boundary  of  Name  of  region  of 


Fig.  3.  The  y-structure  of  the  slice  between  points  P^  and  P,.  in  Fig.  2. 

The  field  "P^Pg"  contains  the  definition  of  the  line  segment 

connecting  the  two  points  P^  and  P^ .  Given  an  x-value,  this  entry 

allows  evaluation  of  the  corresponding  y-value  in  time  0(1). 

Notice  the  entry  R^(!)  under  "Name  of  region."  It  refers  to  a 

tentative  region  which  was  created  at  transition  P^.  In  a  left-to- 

right  scan  of  the  plane  it  will  only  become  known  at  transition  P,. 

that  region  R^  is  to  be  identified  with  region  R^.  The  field 

"Pointer  into  periphery"  will  be  explained  when  we  discuss  the 
r-structure . 


The  r- structure  R 


The  r-structure  is  the  key  to  our  algorithm.  It  integrates  information 
about  the  regions  and  their  peripheries  as  it  is  accunulated  during  the 
unidirectional  sweep.  It  is  initialized  to  be  empty  and  terminates  empty. 

For  any  given  cross  section  it  contains  information  about  exactly  those 
regions  that  are  cut  by  this  cross  section.  Specifically,  with  each 
segment  in  the  cross  section  it  associates  pointers  to  the  two  cyclic  lists 
of  vertices  on  the  boundaries  of  the  regions  which  are  respectively  above 
and  below  that  segment;  equivalently,  with  each  interval  determined  by  two 
adjacent  segments  it  associates  the  cyclic  list  of  the  boundary  vertices  of 
that  part  of  the  region  which  lies  to  the  left  of  the  cross  section.  Figure  4 
shows  the  r-structure  relative  to  a  cross  section  in  the  slice  between  P^ 
and  of  Fig.  2. 


Fig.  4. 


The  r-structure  of  a  cross-section  in  slice  P^P^  of  Fig.  2 
shown  attached  to  the  y-structure. 
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The  description  above  is  a  static  picture  of  R.  Section  2.3  presents 
the  dynamic  picture,  how  R  is  updated  at  each  transition,  and  how  the 
region  boundaries  can  be  printed  to  produce  the  output. 

2.3  The  Region-Finding  Algorithm 

The  algorithm  that  sweeps  the  plane  and  outputs  the  region  boundaries 
has  the  following  simple  overall  structure  (here  X,  Y,  and  R  are  the  three 
data  structures  previously  described): 

Procedure  SWEEP 

begin  X  *■  n  given  points  (sorted  by  x-coordinate) 

Y  *-  (-*,+*) ,  name  of  region  Rq 
R  -  0 

while  X  #  0  do 
begin  P  -  MIN(X) 

TRANSITION (P) 

end 

end 

Procedure  TRANSITION  is  the  advancing  mechanism  of  SWEEP,  and  denotes  all 
the  work  involved  in  processing  one  point  (P)  and  moving  the  "front"  from 
the  slice  to  the  left  of  this  point  to  the  slice  immediately  to  the  right 
of  it;  in  this  process  it  updates  the  corresponding  data  structures  and 
builds  up  the  result  in  an  output  structure.  TRANSITION  will  be  invoked 
exactly  (n+s)  times  and,  since  it  will  be  shown  that  one  invocation  uses 
O(logn)  time,  an  O((n+s)logn)  time  bound  on  the  performance  of  SWEEP  will 
result.  Figure  5  illustrates  the  situation  when  transition  is  invoked: 
for  any  given  segment  s,  A(s)  and  B(s)  are  respectively  the  cyclic  lists 
of  vertices  of  the  regions  bordering  s  above  and  below,  with  the  shown 


orientations . 
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Fig.  5.  The  situation  faced  by  TRANSITION.  The  "next  point"  P  is  of 
one  of  4  types: 

bend  ,  start  ,  end  "Sq  ,  intersection  . 

A(s)  and  B(s)  are  the  cyclic  lists  of  vertices  of  the  regions 
bordering  s  above  and  below,  respectively.  Associated  with  s 
there  are  pointers  to  the  head  of  A(s)  and  the  tail  of  B(s). 

Lists  are  thought  of  as  being  ordered  from  left  to  right;  thus,  for  a 
point  P  and  a  list  L,  "P*L"  denotes  that  P  has  been  added  to  L  as  its  new 
head,  whereas  "L*P"  denotes  that  P  is  the  new  tail. 

Similarly,  for  two  lists  and  Lj*^  denotes  their  concatenation. 
We  shall  now  give  a  detailed  description  of  TRANSITION.  The  function 
INTERSECT (s^jS^)  checks  in  time  0(1)  whether  two  segments  s^  and 
Intersect,  and  if  so,  inserts  the  intersection  point  into  X  in  time 
0 ( log | X | ) . 


procedure  TRANSITION  (P) 
begins  FIND(y[P]) 

case  P  is  "bend": 


(*see  figure  6a*) 


begin  t  segment  beginning  at  P 
h  SUCC(s) 

J i  -  PRED(s) 

if  (f  ^  then  INTERSECT (£ , t) 
if  (h  4  -H»)  then  INTERSECT (t,h) 
A(s)  -  A(s)*P 
B(s)  -  P*B(s) 

Replace  s  with  t 

end 

P  is  "end":  (*see  figure  6b*) 


case  P  is  "end":  (*see  figure  6b*) 

begin  h  SUCC(s) 
t  -  PRED(s) 
i  -  PRED(t) 

if  (f  +  -•)  and  (h  4  +»)  then  INTERSECT (£ ,h) 

Link  A(t)*P*B(s)  and  output  it 
Link  A(s)*P*B(t) 

DELETE (s) 

DELETE (t) 

end 

case  P  is  "start":  (*see  figure  6c*) 

begin  t  -  SUCC(s) 

(l, h)  segments  starting  at  p 
if  (s  4  -®)  then  INTERS ECT(s,X) 
if  (t  4  +•)  then  INTERSECT (h, t) 

INSERT (s ;l) 

INSERTS  ;h) 

Link  B(i)*P*A(h) 

Link  B(h)*P*A(f)  (*B(i) ,A(i) ,B(h) ,A(h)  consist  of  P  alone*) 

end 

case  P  is  "intersection":  (*see  figure  6d*) 
begin  h  *-  SUCC(s) 
t  -  PRED(s) 
l  -  PRED(t) 

if  (X  4  -•)  then  INTERSECT (f , s ) 
if  (h  4  +»)  then  INTERSECT (t,h) 

Link  A(t)*P*B(s)  and  output  it 
A(s)  -  A ( s ) *B 
B (t)  -  P*B (t) 
interchange  s  and  t 

Link  B(t)*P*A(s)  (*A(s)  and  B(t)  consist  of  P  alone*) 


case  P  is  "intersection": 
begin  h  ♦*  SUCC(s 
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Fig.  6.  Illustration  of  the  four  cases  treated  by  procedure  TRANSITION. 

As  a  general  comment  to  the  outlined  algorithm,  all  four  cases  have 
an  analogous  format.  After  segment  s  has  been  located  (Step  1)  in  Y  in 
time  0(log|Y|),  the  other  segments  which  are  candidates  for  intersections 
are  found  in  time  0(1)  (e.g..  Steps  3,  4,  5);  next,  intersections  are 
sought  and,  if  necessary,  X  is  updated  in  time  0(log|x|)  (e.g.,  Steps  6,  7); 
finally,  R  (e.g..  Steps  8,  9)  and  Y  (e.g.,  Step  10)  are  updated,  both  in 
time  0(1)  . 

We  have  seen  in  section  2.2  that  |x|  ■  0(n)  and  | Y |  ■  0(s)  *  0(n  ), 
and  thus  all  the  operations  in  TRANSITION  require  time  0(logn).  Since  the 
algorithm  that  sweeps  the  plane  makes  n+s  transitions,  it  requires  time 
0((rrt-8)logn) ,  thus  confirming  a  prior  claim. 


3.  INTERSECTION  OF  CONVEX  MAPS 


A  map  is  a  planar  graph  G  embedded  in  the  plane:  each  vertex  of  G 
is  represented  as  a  point  and  each  edge  as  a  straight  line  segment  in 
such  a  way  that  edges  intersect  only  at  common  vertices.  A  map  subdivides 
the  plane  into  r  simply  connected  internal  regions  and  one 

external  unbounded  region  Rq.  A  map  is  convex  if  each  internal  region  is 
convex  and  the  complement  of  the  external  region  is  convex. 

Given  two  convex  maps  G^  and  G2  with  n^  and  n^  vertices,  respectively, 
and  s  intersections  of  edges  of  G^  with  edges  of  G^;  we  will  show  that  the 
set  of  these  s  intersections  can  be  computed  in  time  O(nlogrri-s)  and 
space  0 (n) ,  where  n  »  n.+iu.  A  straightforward  application  of  the  plane- 
sweep  algorithm  described  in  [2],  which  does  not  take  advantage  of 
convexity,  would  yield  this  result  in  time  O((rri-s)logn)  and  space  0(n+s). 
The  problem  of  computing  intersections  of  rectangles  with  parallel  edges 
has  been  studied  by  Bentley  and  Wood  [3]  and  McCreight  [5].  Their 
results  are  not  directly  comparable  to  those  of  this  section,  since 
rectangles  may  intersect  even  if  their  edges  don't. 

Let  G^  and  G2  be  two  convex  maps.  For  i  *  1,2,  let  e^^  be  an  edge 
of  G^,  and  assume  that  e^  and  e2  intersect  in  a  single  point  u 
(degenerate  cases  can  be  handled  with  no  difficulty) .  Define  r^  as  the 
region  of  G^  such  that  r^  0  lies  entirely  to  the  left  of  u  (Fig.  7). 

We  now  define  a  plane  domain  U  as  follows.  Let  e ^  j4  be  an  edge  -  if 
it  exists  -  on  the  boundary  of  r^  which  intersects  end  let  e^  be 
analogously  defined.  If  e^  exists,  define  as  the  convex  hull  of  the 
extremes  of  e^  and  e|,  otherwise  is  the  unbounded  half  plane-strip 
orthogonal  to  on  the  side  of  r^;  then  let  U  *  Cl  Uj  (Fig.  7). 
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Fig.  7.  The  domain  U  is  shown  cross-hatched. 

We  claim  that  U  contains  in  its  interior  no  edge,  nor  portion  of  edge, 
either  of  or  of  G2*  It  suffices  to  show  that  Uj  contains  in  its  interior 
no  edge,  nor  portion  of  edge,  of  G^ •  This  is  obvious  when  is  unbounded, 
because  in  this  case  is  contained  in  the  unbounded  exterior  region  of  G^ ; 
when  is  bounded,  then,  by  convexity,  c  ,  and  obviously  the  claim 
holds . 

Consider  now  a  cross-section  Y  at  the  abscissa  of  a  vertex  v  of,  say, 

G^  (Fig.  8a)  and  suppose  that  edges  e^  and  e2»  of  G^  and  G2  respectively, 
are  adjacent  in  Y  and  intersect  at  point  u.  By  the  preceding  argument, 
the  wedge  comprised  between  the  vertical  through  v,  e^,  and  e2  does  not 
contain  edges,  nor  portion  of  edges,  of  G^  and  G2;  hence  we  advance  the 
cross  section  to  include  point  u  to  its  left  (see  Fig.  8b).  This  amounts 
to  exchanging  the  order  of  e^  and  e2  in  Y,  so  that  we  may  proceed  with  the 


(a)  (b) 


Fig.  8.  The  current  cross  section  advances  by  one  point. 

verification  of  whether  e^  and  have  further  intersections  with  edges  of 
and  G^,  respectively.  Specifically,  for  each  intersection  found  two 
new  adjacencies  arise  which  are  potentially  intersection-producing;  each 
such  adjacency  is  then  placed  into  a  queue  for  later  processing.  The 
plane-sweeping  algorithm  will  scan  from  left  to  right  the  vertices  of  both 
maps  and  perform  the  corresponding  deletions  and  intersection  of  edges 
(Steps  5-7  and  9-13,  respectively,  in  the  following  procedure  SWEEP);  for 
any  new  vertex,  at  most  two  adjacencies  are  created  which  may  potentially 
yield  Intersections:  specifically,  if  the  edges  issuing  to  the  right  of 
a  vertex  are  ordered  counterclockwise,  these  adjacencies  may  involve  the 
first  and  the  last  elements  of  this  sequence  of  edges.  After  this 
initialization  of  the  queue  (Steps  12  and  15  of  SWEEP),  for  each  inter¬ 
section  found  the  cross  section  is  dynamically  updated,  i.e.,  an  exchange 
of  order  of  two  segments  takes  places  (Step  23)  and  two  new  adjacencies  are 
examined  for  possible  addition  to  the  queue  (Steps  20-22).  At  any  point 
in  the  execution  of  the  algorithm  the  cross  section  corresponds  to  a  curve 
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in  the  plane  (not  at  straight  line,  necessarily!)  which  bounds  the 
discovered  Intersections  to  their  right  (Fig.  9). 


Fig.  9.  Structure  of  cross  sections  in  the  proposed  technique.  Edges 
of  are  shown  in  solid  lines,  those  of  in  broken  lines. 

A  description  of  the  procedure  is  given  below.  For  each  vertex  v  (either 
of  G^  or  G2)  we  denote  as  L(v)  and  R(v)  the  sets  of  edges  incident  to  v 
and  lying,  respectively,  to  its  left  and  to  its  right.  The  y-structure 
Y  is  a  dictionary,  while  the  x-structure  X  is  simply  an  array;  I  denotes 
the  set  of  the  intersections  found;  Q  is  a  queue  (implemented  as  a 
linear  list),  for  which  '**  Q"  and  "Q  «"  denote  the  operations  "remove" 
and  "add."  For  simplicity,  the  algorithm  omits  consideration  of  the 


r-structure. 
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procedure  SWEEP 


1 .  begin 

for  each  vertex  v  of  G^UG„  do  sort  counterclockwise  the 
terms  of  L(v)  and  R(v) 

2. 

sort  the  vertices  of  G.UG?  by  increasing  abscissa  and  place 
them  in  X[l:n] 

3. 

Y  -  0,  I-«,  Q  -  0 

4. 

for  i  *-  1  until  n-1  do 

5. 

begin  while  (L(XfiD  4  0)  do 

6. 

begin  e  •”  next  edge  in  L(Xfil) 

7. 

DELETE (e) 

end  (*all  edges  incident  to  X[i]  from  the  left 
are  deleted*) 

8. 

s  -  FIND(y(X[i ]) ) 

9. 

while  (R(Xfil)  4  0)  do 

10. 

begin  e  next  edge  in  R(Xfil) 

11. 

INSERT  (s;e) 

12. 

if  (s  and  e  belong  to  different  maps)  then 

Q  «  (s,e) 

13. 

s  •-  e 

end  (*all  edges  issuing  from  X[i]  to  the  right 
are  inserted*) 

14. 

s'  -  SUCC(s) 

if  (s  and  s'  belong  to  different  maps)  then  Q  c  (s,s 
(*at  most  two  pairs  of  edges  have  joined  Q*) 

15. 

16. 

while  0  4  0  do 

17. 

begin  (e^e,,)  «  Q 

18. 

if  (e^  and  e^  intersect)  then 

19. 

begin  I  *-  lU(e^,e„) 

20. 

e'-  PRED(e^)  ,e"  -  SUCC(e2) 

21. 

if  (e'  and  e2  belong  to  different 
maps)  then  0  *  (e',e„) 

22. 

if  (e^  and  e"  belong  to  different 
maps)  then  Q  (e^,e") 

23. 

exchange  (e.,e„)  in  Y 

end 


end 


end 
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Notice  that  intersections  are  discovered  only  when  marching  on  an 
edge  from  lef t-to-right;  where  an  edge,  is  deleted  (Step  7)  all  inter¬ 
sections  of  this  edge  have  already  been  found.  Moreover,  when  edges  are 
deleted  no  intersection-producing  adjacencies  are  generated;  indeed  if 
e  (being  deleted)  is  incident  to  v  from  the  left,  either  there  is  some  e' 
incident  to  v  from  the  right  which  maintains  separation,  or  v  is  the 
rightmost  vertex  of  its  map:  but  in  this  case,  no  adjacencies  between  edges 
of  and  G^  are  possible  to  the  right  of  v.  The  running  time  of  the 
algorithm  is  easily  shown  to  be  O(nlogn).  Indeed  such  is  the  running  time 
of  Steps  1  and  2;  loop  5-7  uses  O(logn)  time  for  each  vertex,  and  is 
executed  n  times;  similarly  does  loop  9-13.  Notice  that  Steps  12,  15, 

21,  and  22  each  use  0(1)  time,  and  they  are  collectively  executed  0(s) 
times;  finally,  loop  16-23  involves  0(s)  manipulations  of  the  y-structure 
(Steps  20  and  23)  each  of  which  uses  0(1)  time,  since  pointers  to  e^  and 
are  available.  Thus  we  conclude  that  our  map  intersection  algorithm  runs 
in  time  0(nlogn+s).  The  0(n)  space  bound  is  obvious. 
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4.  VARIANTS  OF  THE  PLANE -SWEEP  ALGORITHM 

In  order  to  assess  the  generality  of  the  plane-sweep  algorithm,  we 
cast  the  two  instances  used  in  sections  2  and  3  into  a  common  frame. 

Both  algorithms  have  the  following  structure: 

Algorithm  SWEEP: 

1.  Initialize  x-structure 

y- structure 
(perhaps  other  data) 

2.  while  x-structure  not  empty  do  TRANSITION 
where  TRANSITION  is  of  the  form 

1.  P  *  (x,y)  *-  remove  next  point  from  x-structure 

2.  with  y  locate  an  interval  in  the  y-structure 

3.  compute  i^  new  intersections  and  process  these. 

The  y-structure  is  identical  for  both  algorithms.  It  stores  the  current 
cross-section  consisting  of  0(n)  entries  in  a  data  structure  that  supports 
the  operation  FIND  in  logarithmic  time,  and  the  operations  PRED  and  SUCC 
in  constant  time. 

The  x-structure  is  rather  different  for  the  two  problems  we  have 
discussed.  The  simple  case  is  illustrated  by  the  convex  map  problem: 
all  relevant  transitions  are  known  a  priori:  the  n  =  nj+nj  vertices  of 
the  two  given  graphs.  After  they  have  been  sorted  they  can  be  stored 
in  any  static  data  structure  suitable  for  sequential  processing  (i.e., 
the  operation  NEXT  takes  constant  time)  -  for  example  an  array.  The 
reason  is  that  the  intersections  being  computed  can  be  processed 
entirely  (an  0(1)  operation)  when  they  are  encountered;  hence  they 
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don't  need  to  be  considered  to  be  transitions,  i.e.,  they  don't  need  to  be 
stored  for  deferred  processing  and  retrieved  from  the  x-structure 
(an  O(logn)  operation).  By  contrast,  in  the  regions-of-a-polygon 
problem,  a  computed  intersection  must  be  treated  as  a  transition,  to  be 
stored  and  retrieved  from  the  x-structure.  This  requires  a  dynamic  data 
structure,  which  supports  the  operations  MIN  and  INSERT,  and  cannot  be  as 
efficient  as  a  static  data  structure.  The  mere  fact  that  operations  on 
the  x-structure  now  require  logarithmic  as  opposed  to  constant  time, 
however,  would  not  affect  the  asymptotic  time  requirement  of  the  algorithm, 
since  this  access  time  gets  absorbed  in  the  nlogn  term.  The  difference 
between  O(nlogrH-s)  and  O((n+s)logn)  is  merely  due  to  the  fact  that  n+s 
transitions  move  through  the  x-structure  as  opposed  to  n. 

This  comparison  is  summarized  in  the  following  table: 


regions -in-a- loop 

convex  maps 

SWEEP: 

1.  Initialize 

sort:  0 (nlogn) 

sort:  0 (nlogn) 

2.  while -do  TRANSITION 

(n+s)  times 

n  times 

TRANSITION: 

1.  P  remove  from  x-structure 

MIN:  0 (logn) 

NEXT:  0(1) 

2.  locate  interval  in  cross 
sec tion 

FIND:  O(logn) 

FIND:  0(logn) 

3.  process  i^  new  intersections 

i  =  0,1  or  2 

i  times  0(1) , 

p 

INSERT:  O(logn) 

where  2  i  =  s . 

P 

TOTAL: 

0 (nlogn) 

0 (nlogn) 

+  (nt-s)O(logn) 

+  s  *0  (1) 

The  two  algorithms  presented  can  easily  be  combined  to  compute  the  regions 
of  the  intersection  of  two  arbitrary  maps  (non-convex)  in  time  0((n+s) logn) . 
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In  order  to  do  this,  however,  the  classification  of  points  into  the  4 
categories  (sect^on  2.3)  must  be  changed  to 

deal  with  one  general  type  of  point  where  an  arbitrary  number  of  edges 
meet.  This  modification  resolves  the  problem  of  degeneracy  mentioned  in 
section  2.1.  An  intersection  between  more  than  2  edges  in  the  same  point 
is  simply  treated  as  a  vertex  of  high  degree.  By  means  of  the  same 
generalization,  the  regions  of  the  intersection  of  two  convex  maps  can 
be  computed  in  time  O(nlogrrfs). 
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